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Spectra of one-particle and collective excitations in narrow-band ferromagnets with unquenched 
orbital moments are calculated in various theoretical models. The interaction of spin and orbital 
excitations with conduction electrons results in the damping of the former which, however, turns 
out to be rather small; therefore, apart from usual spin waves, well-defined orbitons can exist. 
Non-quasiparticle states occur in the electron energy spectrum near the Fermi energy due to this 
interaction. The criteria of stability of the saturated spin and orbital ferromagnetic ordering are 
considered. Possible effects of orbital ordering in magnetite and in colossal magnetoresistance man- 
ganites are discussed. 
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I. INTRODUCTION 

Orbital ordering in strongly correlated systems is a subject of interest for a rather long period pj , this interest grow- 
ing considerably last time (see, e.g., recent papers [2, 0, La M ) ■ F° r orbitally ordered magnets the energy spectrum 
is more rich than for purely spin- moment magnetic materials since it contains additional excitation branches ("or- 
bitons"). Relevance of the orbital degrees of freedom for physics of the colossal magnetoresistance (CMR) manganites 
0,0, HSU is one of the reasons for the revival of this field. A large orbital contribution to the total magnetic moment 
has been observed recently by the X-ray magnetic circular dichroism (XMCD) method for magnetite Fe304 [l0|. The 
latter substance (the most widespread natural iron compound) is probably one of the most intriguing magnetic ox- 
ides because of controversial experimental information concerning its ground state and so called "Verwey transition" 
[HD3 (see also a discussion in 

Refill). 

The magnetite is also a promising material: it is a half-metallic ferromagnet 
(more exactly, ferrimagnet) with a high Curie temperature and perspectives of applications in spintronics |15| . Since 
current theoretical description of the half-metallic ferromagnetism (for a review, see Refs lla. Il6|) ignores the orbital 
degrees of freedom, it should be generalized to include them and to take into account probable orbital ferromagnetism 
in the magnetite and, possibly, in other strongly correlated conducting ferromagnets. Such a generalization is the aim 
of the present work. 

The paper is organized as follows. In Section 2 we formulate the models of a narrow-band ferromagnet including both 
spin and orbital degrees of freedom and discuss their relevance for the colossal magnetoresistance (CMR) manganites 
and for the magnetite. In Section 3 the Bose-type Green's functions are calculated which determine the spectrum and 
damping of spin waves (magnons) and orbitons. In Section 4 the Fermi-type one-particle Green's functions are derived 
and peculiarities of the electron energy spectrum owing to the electron-magnon and electron-orbiton interactions are 
considered. Using the expressions obtained the conditions for stability of saturated spin and orbital ferromagnetic 
ordering are investigated. In Section 5 a general picture of the half-metallic ferromagnetism in the presence of the 
orbital degrees of freedom is discussed. 

II. MODELS OF A CONDUCTING FERROMAGNET WITH ORBITAL DEGREES OF FREEDOM 

Previous works concerning the role of the orbital degrees of freedom in strongly correlated magnetic systems treated 
mainly the insulator case, in particular, orbital ferro mag netism accompanied with spin antiferromagnetism [J], orbital 
antiferromagnetism H, or orbital liquid 0- In Rcf 11 a mean-field phase diagram of an insulator including orbital 
and spin ferromagnetism region was obtained. In contrast with these works, we consider here conducting ferromagnets 
with both spin and orbital degrees of freedom. 

We start with the many-electron system with two ground terms of the d n and d n+1 configurations, T n = {SL} and 
r n+ i = {S' L'}. The corresponding periodic Anderson model that describes the electron hopping with the transitions 
between two these states only reads 

H = H ° + X! ^ C L C k,m + X { V ^<ia m ^ a ^ra + h.C.) (1) 
kcrm kimcr 

where TLq is the Hamiltonian for strongly correlated d-electrons described by the operators au am , ikm is the band 
energy for conduction electrons, Vk is the matrix element of hybridization. It is suitable to use the representation of 
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the Fermi operators for the localized electrons in terms of the many-body atomic quantum numbers |17| 

4™ = (» + !) V2 E G 1 L '<KS.^ (S'L W, SLM») (2) 

where Gf/' are the fractional parentage coefficients, Xi(a(3) — \ia)(i/3\ are the many-electron Hubbard operators 
[l8j which satisfy the relations at a site i 

X t (a{3)X t ( 7 e) = S^X^ae), £ X^aa) = 1. (3) 

a 

By a canonical Schricffcr- Wolff transformation the model Q is reduced to the broad-band s — d exchange model 
with the exchange parameter I sd ~ V 2 /A (A is the difference of the term energies). In a general case, we have the 
s — d exchange model with two magnetic configurations. This model was used to investigate the electron spectrum 
and calculate the Kondo temperature in RefslrH. IT^ . Here we introduce a further simplification assuming that only 
one of the competing configurations has non-zero orbital moment L — I. This assumption holds for the magnetite 
which has d 5 and d e ground-state configurations for Fe 3+ and Fe 2+ respectively, the first configuration having zero 
orbital moment. A similar situation takes place for the CMR manganites (with d 3 and d 4 configuration for Mn 4+ and 
Mn 3+ ) : due to the relevance of ti g — e g crystal-field splitting the former configuration corresponds to the completely 
filled t2 g band with L = 0. 

According to RefsIT^. Il9l the s — d exchange Hamiltonian takes the form 

n sd = -I sd £ Cl'^Cl'^ ^XiiSLMm, SLM"m)ci mCia , m . (4) 
Rearranging the product of the Clebsch-Gordan coefficients C with the use of 6j-symbols, 

E n s'M' r<s'M' _ 3\/2 / 25^ + 1 \ 1/2 J 1/2 1/2 1 1 n SM n h° fz\ 
°SM,Ia°SM",ia' ~ ~ \ 2S + 1 \ S S' S' f 2^ ° ' SM " ^ 'W ( ) 

M> \ ' / V J q= _ 1 

and taking into account that C§^,, lq are the matrix elements of spin operators, we obtain the s — d exchange 
Hamiltonian with orbital degrees of freedom in the form 

+ n d (6) 

ktrm iaa'ra 

where 7i d is the Heisenberg Hamiltonian of the localized-spin system, er are the Pauli matrices, the s — d exchange 
parameter / is renormalized in comparison with I sd in Eq.Q). A similar one- impurity model was discussed for the 
Mn 2+ ion (orbital singlet problem) |2(J, [21j . 

In contrast with the case of rare-earth elements where the orbital moments are not quenched, but act as a part 
of the total orbital moment (the Russel-Saunders coupling [171 Il9l 22]), for 3c?-electron systems we have to consider 
the case of strong crystal field where the excitation spectrum is characterized by spin and orbital quantum numbers 
separately. However, in the case under consideration the quantum number m enters the model in a very simple way 
(the Hamiltonian is diagonal in the orbital indices) and can label both spherical and cubic harmonics. Of course, a 
relatively simple form of the Hamiltonian © results from our assumption that only one configuration has a non-zero 
orbital moment. Fortunately, it is the case that is most interesting for real systems, as it was mentioned above. 
Thus a complicated mathematics involving the orbital Clebsch-Gordan coefficients is relevant only when we have two 
many-electron terms with L =/= 0. 

Note that in the case of j — j coupling (e.g., 5/-electron systems) where the electron states \jm) are characterized 
by the projection of the momentum j, we have the Hamiltonian 

km i 

3i = ^jm,la C L C ™'' CO 

mm' 

Provided that the crystal field Hamiltonian is diagonal in m, the model © does not mix (at least in the lowest 
orders of perturbation theory) the states with different m, see Section 4. Therefore we focus in the present work 
mainly on a more rich and complicated narrow-band case (which is also most interesting for real transition metal 
compounds) . This case should be described by a two-configuration Hubbard model where both conduction electrons 
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and local moments belong to the same (i-band, the states with n+1 electrons playing the role of current-carrier states. 
After performing the procedure of mapping onto the corresponding state space, the one-electron Fermi operators for 
the strongly correlated states a\ lljm are replaced by many-electron operators according to Eq. J2J) . Taking into account 
the values of the Clebsh-Gordan coefficients which correspond to the coupling of momenta 5 and 1/2 we obtain 

kcrm 

Here we have redefined the band energy by including the many-electron renormalization factor, 

t km {n + l)(Gfl°) 2 /(2Z + 1) -» < km , (9) 



and 



s 



M=—S 
S 



9tm= E y^sTT X ^ S ~ 1 ^ M + Z-SMm), S' = S-l/2, 



^ v 5 + 2s + i l xt{s + 1/2 ' M + V SMm) > s ' = s + 1 / 2 ( 10 ) 

M=-S 

where \SMm) are the empty states with the orbital index m, \S'M') are the singly-occupied states with the total 
on-site spin S' — 5 ± 1/2 and its projection M', a = ±. 

We see that the two-configuration Hamiltonian is a generalization of the narrow-band s — d exchange model with 
|J| — > oo fT?l l23l]. In the case where the configuration d n+1 has larger spin than the configuration d n , we have the 
effective "s — d exchange model" with / > 0, and in the opposite case with / < 0; it is worthwhile to remind that our 
band energy is renormalized according to © in comparison with the s — d exchange Hamiltonian (cf. the factor of 
1/2 for the narrow-band Hubbard model without orbital degeneracy |23j). 

In the case of j — j coupling we have 



9i(7m * 9im 



lm — 2l CjM <jm Xi(j'M' , JM) 



MM' 



The specification of the Hamiltonian JHJ) for the case of a saturated ferromagnetic state (where only maximum spin 
projections give a contribution) reads for I < (5' = S — 1/2) 



h — ^2 <k« 

km 



25 v a -+ m v+m,a . 1 y»- m -y — m,a , y 25(2S 1) 



X k A k + OC i -I A -k A k H nr, | 1 A -k A k + n - C - 



25 + 1 - k k 25 + 1 ~ k k 25 + 1 



(ii) 



Here | + m) and | — m) are the on-site states without excess electrons and with localized spin projection 5 and 5—1, 
respectively; \a) = \S — 1/2,5— 1/2) and \b) = \S — 1/2, 5 — 3/2) are the states with one excess electron and the total 
spin projection 5 — 1/2 and 5 — 3/2, respectively; X?r are the Fourier transforms of the Hubbard operators Xj,(a, (3). 
For the case where / > (5' = 5 + 1/2) we have 



H = ^2 
km 



2"? 1 9^ 

x - k * k + 25TT x - k * k + 25TT x - k Xk + V^5+T (x - k * k +h - c - } 



(12) 

where \u) = \S + 1/2, 5+ 1/2) and |y) = |5 + 1/2, 5-1/2) are the singly-occupied states. 

We will consider also the simplest model where one of the configuration corresponding to the "current carriers" has 
zero spin and orbital moments (non-magnetic holes): 



H=J2 t±mXX m Xr'°- (13) 

kfrrn 

This is a formal generalization of the standard narrow-band Hubbard model in the X-operator representation to 
include the orbital degeneracy. Apart from the m-dependence of the band energy £km, the Hamiltonian l|13fl includes 
spin and orbital degrees of freedom on equal footings. However, for finite Hubbard repulsion U the symmetry of 
spin and charge channels is in fact lost owing to the Hund interaction, so that the state with orbital ferromagnetism 
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and spin antiferromagnetism can become favorable in comparison with the spin antiferromagnetic state [l| . Note the 
equivalence of the model (|13|) with the replacement t km — ► ikm/2 and the model (|llf> for S = 1/2. 

For orbitally degenerate case the subbands t km are connected by the point-group symmetry transformations and 
have identical densities of states. However, in real situations usually the Jahn- Teller lattice instability takes place 
which destroys the exact orbital degeneracy Due to the Jahn- Teller distortions, both additional on-site crystal- 
field splitting arises (which means the shift of the centers of the subbands) and hopping integrals are changed (which 
means different k-dependences). Thus we will assume further that all the functions tkm are in general different. 



III. COLLECTIVE SPIN AND ORBITAL EXCITATIONS 



The spectrum of electron and spin excitations in the above-discussed narrow band models can be calculated similar 
to the papers [23,12^,12^1, the difference being in the occurrence of transitions into orbital states. 

Even in the narrow band case where small interaction parameter is absent we have an hierarchy of energy scales 
owing to both formal expansion parameters (inverse nearest-neighbor number 1/z, quasiclassical parameter 1/S) and 
small current-carrier concentration. Therefore the width of the "magnon" band is much smaller than characteristic 
electron energies. A similar approach was developed earlier for degenerate ferromagnetic semiconductors pflj . 

First we calculate the retarded commutator Green's function describing spin and orbital excitations in the simplest 
model (JT3J, 

<WM = ({X™«r' m '\X^'' am )) m Imcj > 0. 

To derive the equation of motion, we use the commutation rules for X -operators, which follow from Eq.Q. We 
obtain 

vG& m ,(w) = N am - AW + - ^{{X^™' 'xr'°|r'; >m )) u . (14) 

k 

In the next equation of motion we make a simple decoupling described in Refl25l which corresponds to the first-order 
term in 1/z-perturbation theory (each order in 1/z corresponds formally to a summation over an additional wave 
vector) : 



[V — ^kam + fi k-q<T'm'jU A q ^k A k l A -q //w ~ n kam ~ n k-qcr'm' 

+ (ikm^kcrr?x ~ tk-qm'Bk-qir'm')'j'qmm'( w )' 

Here the "Hubbard-I" energies E kcrm and the occupation numbers nko-m are given by the "Hubbard-I" |28| expressions 

Ekam = tkm(N am + No), N a — (X™ a ) , 
),crm V (T 

k 

f(E) being the Fermi function. Then we have 

(^k-qm' ^Icm)(%f7m ~~ ^-k— qc'm') 



n kan = ( X °jrx° m '°) = (N am +N )f(E ]i<7m ), (16) 



(^k-qm' ^km) (^km^ktrm ^k-qm' ^k— qcr'm' 

k 



^ ^ U) ^kcrrri ^k— qcr'm 



(17) 



In the saturated ferromagnetic state (N am = 8 a +N m ) and for small hole concentrations c = N the pole of the Green's 
function 1)17(1 yields a simple result for the spin- wave frequency 



LO+- , 
qmm 



y^ftk-qm' — tk m )f(tk m ). (18) 



In the case of saturated orbital ordering (N am — 8 m+ N a ) the orbital excitation frequency w^ m , (m = + ^= m') is 
given by the same expression i[18jl. Provided that the states m' and m are split by the crystal field, there is a gap in 
the corresponding orbiton spectrum branch. Assuming for simplicity that the band bottom corresponds to T-point 
k = 0, we derive 

<W - c(W " W) + cA, A = - W (19) 
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It should be noted that the gap c^o = is proportional to the current carrier concentration, but not just equal to 
the crystal-field splitting A, in contrast with the case of weak crystal field typical for the rare-earth systems [2^| . 
Actually, in our case characteristic orbiton energies for finite wave vectors have the same order of magnitude as the 
magnon energies. 

The magnon damping can be obtained similar to Refl23l bv calculating higher-order terms in 1/z, 

7qmm'(w) = 7T ^ (tk-qm' — *km) 2 ["k-q+pm" (1 — "km) + Nb {<Jj p ) (rik-q+pm" ~ "km)] S(uJ + £ km — ^k-q+pm" — ) 
kpm" 

= 71 X! fa-qm' - t km ) 2 (uj p - UJ) 9 ^" 1 [N B {UJ P ) ~ N B (u p - u)]6(t km - tk-q+pm" ) (20) 

, ^km 
kpm" 

with Nb(w p ) being the Bose function. The damping turns out to be finite at T = 0, unlike the case of the Heisenberg 
ferromagnet. However, it is proportional to the frequency ui and contains formal small parameters, which leads to 
very small value |23| . Thus both spin-flip and non-spin-flip orbital excitations are well defined. This fact is non-trivial 
in a narrow-band case (a similar situation takes place also in the antiferromagnetic case, cf. Ref@). 

In the model (JHJ the results for the "magnon" spectrum are more complicated and depend in a non-trivial way on 
the spin value. The calculations for the ferromagnetically saturated state can be performed similar to Refl23l For the 
case S 1 = S — 1/2 we write down the sequence of equations of motion 

-G+ mm ,M = N m + £[(t k _ qm , - 2St km )((X^'x+ m ' a \XI™'> 

k 

+ ^25(25 - l)ty m ,{{X^+ m 'x+l l *\XZ™'- +m ))u]- (21) 

The equation for the first Green's function in the right-hand side is obtained as above (see Eo.l|15|l). and for the 
second one we have 

(w - ti m >N m )((X a J+ m 'x+^ b \XZ^'+ m )) u = ^2S(2S T)(X a J+ m X+ m ' a )t Um ,G+- m ,(u;) (22) 



where 

25 



Km 2S+1 

Then the excitation energy in the leading order in 1/z reads 



tkm- (23) 



<mw = 4 £[(*k- q w - 2Stljm m ) + (2-5 - l)C'/(*km')] ( 24 ) 



qmm' 2^ 



(more strict calculations in the absence of orbital degeneracy |23| yield a closed integral equation for the magnon 
spectrum) . At the same time, the frequency of the non-spin- flip orbital transitions is obtained similar to Eq. i|18|) and 
is given for a saturated orbital-ordered state by the expression 



qmm 



' — y^ftk-qm' *km)/(*km)) (25) 



k 



so that a dependence on S is absent. 

For S' = S + 1/2 we have the equations of motion 



(w - t km N m )((X^ X+ m < u \XI™ > +m )) w = -t km (X^ m X+ m -)G+ mn ,(L,), (27) 



tu m N m - t k+ « m ,N m ,/(2S + l)]((A! k +m A k + ;;' "|Jf I™'- = ^^||-(^ k +m A k + ™-")G+ mm ,(c). (28) 
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Then we derive 

k 

The orbital frequency for a saturated state reads 



qmm 



/ — ^(^k+qm' — ^km)/(ikm)- (30) 



It should be mentioned that, instead of the difference tk m N m — tk+qm'-^m', a resolvent occurs in the denominator 
in the right-hand side of Ea. I|29fl at more rigorous calculations which do not use the 1/z-expansion (cf. Ref l23T) . In 
particular, in the case of ferromagnetically saturated spin and orbital ordering (all the sites without excess electrons 
are in the orbital state | + m) = | + +) with some to = +) we obtain the following exact result for the magnon pole 
(m = rn'): 



+- 

qm 



k 



2*5 + (tkm ^k+qm)-Rkm(<- i -'qm ) 
1 — /(^k+pr, 



i? km H = E ' fpggi _ , (32 ) 

^ T" tkm ^k+pm ^pm 



P 



Note that the imaginary part of the expression l|31|l describes the magnon damping. 



IV. ELECTRON SPECTRUM AND INSTABILITIES OF THE SATURATED STATE 

In the broad-band s — d exchange model (JBJ) we can perform a decoupling in spirit of a ladder approximation (cf. 
ReflU to obtain the Green's function for the case of zero temperature T = 

fe™l4™))E = [E- t ktTm - ^(E)}- 1 (33) 
where t^ arn = tkm — crlS is the mean-field (Hartree-Fock) spectrum. The electron self-energy reads 

E -< E > - ttttse? (34) 



«L(£) = E = E g < 35 > 

^ — Ik-qJ.m + w q q ^ — tk-qtm " w q 

with raj = /(tko-m)- Imaginary part of (|34() describes the non-quasiparticle states which are due to non-pole 
contributions (branch cut of the Green's function) 0,01- ^ n the saturated ferromagnetic case n k ^ for / < 
and ral ^ for I > 0, and these states with a — —sign/ lie below (above) the Fermi level, respectively. The 
corresponding densities of states for a twofold orbital-degenerate band are shown in Figs. 1,2. One can see that in 
the model @ the subbands with different to are uncoupled (we do not consider their mi xing owing to off-diagonal 
crystal field), so that the picture of the half-metallic ferromagnetism state discussed in Refs ll5ll24l is not qualitatively 
changed, although the number of excitation branches increases. 

A more complicated situation takes place for the narrow-band limit, a formal reason being in occurrence of inter- 
subband transitions owing to non-trivial commutation relations for many-electron X-operators. First we consider 
the one-particle spectrum in the simplest model l|13() . We treat the case T = and restrict ourselves to the case of 
spin- and orbital-saturated ferromagnetic state (all the singly occupied on-site states have the spin projection a = + 
and the orbital state m = +). Since this ground state is non-degenerate, the one-particle spectrum at small hole 
concentrations c — Nq can be investigated in a strict way similar to Ref(2^ We have to calculate the retarded 
anticommutator Green's function 

G k<Tm (E) = {(Xr'°\X°X m }}E, lm£>0. (36) 
The current carriers with \am) = | + +) propagate as free ones, 

G k++ (E) = (E-t k+ )-\ (37) 
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In the case \am) ^ |++), to extract the operators of spin and orbital excitations, it is convenient to use the kinematical 
relation which follows from Eq.|j3J) 

xr ,o = £ xr , ++x ++,o. (38) 
q 

Then we have 

G^ m {E) = F ^™(E), Fkqo-m (E) = ((X^ m - ++ X++^\X°X m )) E . (39) 



After a natural decoupling procedure (cf. Ref.l24f) we obtain the integral equation for the function F 

(E - i k - q+ + u+l m )F^ am {E) = n k _ q [l - (t k _ q+ - t km ) ^ F kpam (E)} (40) 

p 

where n k = (X°j£ + X^ + '°) = ,/(f k +). Solving this equation we derive 

G k am(E) = L (41) 

E — tkm + 25/ Ri iam (E) 

where 25 = 1 — c is the average spin magnetization, 

q & - t k -q+ + Wq+m 

Thus the Green's function with minority spin (orbital) projection has non-pole structure, the non-quasiparticle states 
below the Fermi level being of crucial importance to satisfy the sum rule 

f dEf(E)g am (E) = T£ {X °X m X™>°) = N , 
J k 

9a m (E) = --ImVG kffro (£), (43) 

k 

irrespective of a, m. 

The spectrum picture is presented in Fig. 3. To simplify numerical calculations, we average the resolvent i? kcrm (-E) 
in 141|) over k, 



Rkam(E) -> i? CTm (^) = J dLjK am {u)Y^ 



" k 



£ - i k+ + w 

(this approximation yields the correct behavior near the Fermi level, cf. Refl2(j, although yields an irrelevant unphysical 
shift of the band bottom by the maximum magnon frequency). We use the semielliptic magnon density of states 
K am {uj) which is proportional (with the corresponding shift) to the bare electron density of states according to 
the band edge uq being determined by the crystal-filed splitting. When neglecting spin (orbital) dynamics, 
Imi? kcrm (i?) has a jump at the Fermi level. For gapless magnons lmR am (E) (this quantity is also shown in Fig. 3) 
vanishes at Ef according to the law (E F - Ef/ 2 \2$. The "Kondo" singularity of ReR am (E) owing to the Fermi 
function n k results in that the dependence of g a m{E) near Ep is considerably more sharp than just of lm.R am {E). 
The contributions of orbitons with a gap spectrum demonstrate a threshold energy equal to the gap, that is, in the 
presence of the orbital splitting lmR am (E) starts from Ep — uiq. At the same time, the threshold energy decreases 
strongly for the total density of states g am {E). Besides that, the orbital splitting results in a height increase of the 
9am(E) peak below the Fermi level. 

The instability of the saturated state (N am ^ for \am) ^ | + +)) with increasing c corresponds to occurrence of 
a non-zero spectral density above Ep since 

' dE[l - f{E)]g am {E) = N am . (44) 
This spectral density comes from the quasiparticle pole determined by the equation 

RkamiEp) = — . (45) 

max r km — Ef 
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The analytical estimations are complicated due to the logarithmic singularity in the resolven t Ry r, m (E). It should be 
noted that a simple cutoff of this singularity at a characteristic magnon frequency (cf. Ref. |22|) is insufficient since 
the energy dependence of the resolvent plays an important role (Imi? ko - m (E 1 ) oc (Ep — E) 3 / 2 and the maximum in 
Rci? ko . TO (_E) does not correspond to Ep). Thus an integration with the magnon spectral density is required. As it was 
demonstrated by the corresponding numerical calculations for the Hubbard model [26j, the critical concentration of 
holes Ccrit makes up about 20-30% for a number of crystal lattices. To calculate the second critical concentrations c' clit 
(a transition from the non-saturated ferromagnetism into the paramagnetic state), more advanced approximations 
taking into account longitudinal spin fluctuations are necessary 26]. Note that in a more realistic case of finite on-site 
interaction the conditions for a ferromagnetically saturated state can be considerably more strict. 

Provided that the spectra i km with different m are connected by crystal symmetry transformations (e.g., states 
corresponding to the irreducible ti g representation in a cubic field), orbital and spin instabilities of the saturated 
state coincide. However, if the orbital subbands are split by the crystal field A (or, in other words, an orbital 
pseudomagnetic field is present), the criteria for the orbital and spin- flip instabilities can become different owing to 
the difference in the spectra i km and i k + and the gap in the spectrum w q + m . 

The analytical calculations of the one-particle excitation spectrum for the model (|IIJ) are performed in a similar 
way. We restrict ourselves again to the case of a small conduction electron concentration c. In the case S' = S — 1/2 
spin-down electrons have the spectrum t k+ (see Ea. (l23|) ). For the function 

F kq (£) = {(X--++X+±£\XX + ))e (46) 
we obtain after some simplifications the integral equation 

(E - ^_ q+ + <+ + )F kq (£) = n k _ q [l - ^(t k _ q+ - t£ + ) J2 F k P (£)] ( 4? ) 



p 



with n k — /(t k+ ). Then we derive 



2S 1 

o m (B) = pt+-|*V>. = E -^ + (1S - c)/Rum < 48 > 

where -R ko . m (-E) differs from (|42|l by the replacement i k+ — » f J , . The condition for the spin-flip instability corresponds 
to occurrence of the pole of the expression l(48j) above the Fermi level. As follows from numerical estimations, such a 
condition can be hardly satisfied for 5 > 1/2 (although the magnon frequency is also more soft in this case because 
of the factor (2S)- 1 in Eq.JSJ). 

The orbital instability is essentially governed by the Hamiltonian of the structure (|13|l (the first term in Eg. 1111(1 ) 
and is determined by the Green's function 

^'^■- B^rR^TCT - (49) 

Since in such a case the factor of 25 at 1/R*(E) is absent, orbital ordering becomes unstable at some electron 
concentrations, and we can treat the state with saturated ferromagnetic spin ordering, but destroyed orbital ordering. 

In the case S' = 5+1/2 (model 112(0 spin-up electrons behave as free ones, and we have for the spin-down Green's 
function 



[E-t lc+ /(2S+l)]((X^ v \X v j+ + )) E = l- c +^-!^£ tk _ q+ L kq (£), (50) 

q 



L kq (£) = ((A++'-+A+_ + q tl |Al k ++ )) f; , (51) 



25 . „,, ++i1) , v v,++\ 



(£-t k _ q+ -u; q )£ kq (£;) - ^^^ k _ q+ ((A™|A! k + - 

-J^tk-p+fl-nk-pJWtf) ( 52 ) 
p 

where n k = (X^ + X+ + ^ u ) - /(* k+ ), Wq — w q ^_|_. Solving this system we obtain the expression which also has a 
typically "non-quasiparticle" form (without quasiparticle poles, at least, at not too large current carrier concentration) 

((x ++,„ |x!> + +)) = 2S+l (53) 

E - t k+ + 2S/R k (E) 
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Mm =J2 e h n " q u ■ (54) 

Note that the result (|53[) is rigorous for the empty conduction band (spin-polaron problem) and is in agreement with 
(|33[I - H35[I at I —* +00. As follows from the spectral representation, n£* = (X^ + X^ + ' v ) = (provided that the 

Green's function l|53|l has no poles) since ImR^E < Ep) — because of the factor 1 — rtk-q- Physically, there is no 
spin-flip instability in the s — d model with large I > since spin-down electrons are not present in the ground state. 
The corresponding spectrum picture is shown in Fig. 4 

Note that the Green's function (|53|l itself does not determine an instability of the saturated state since the spectral 
representation gives 

- I f dE[l - f(E)]ImJ2((K m ' V \X V -r))E = N am . (55) 
To treat the instability we have to consider the Green's function with a = — , 

{{x - m , V]x v,- m))E = J2 {{X -m, ++X ++^ xVj - m)) ^ (56) 

q 

which is proportional to the distribution function n£* and is expressed in terms of the resolvent 

^k* 

-Rk+m = 2^ p 7^; —7— (57) 

q & ~ *k-q + W q 

where i£* is the spectrum of new quasiparticles which can occur with increasing electron concentration. 
The orbital instability is determined by the first term in the Hamiltonian l|12|) . 



= w^r^m^E-y ,58) 



and is again S- independent. 



V. CONCLUSIONS 



In this paper we have presented the picture of excitation spectrum in a saturated conducting (half-metallic) ferro- 
magnet with orbital degrees of freedom within the framework of simple many-electron models. It is proven that the 
half-metallic ferromagnetic state does exist in the presence of orbital degeneracy. We have focused our consideration 
on strongly correlated systems with pronounced term (multiplet) effects. In such a situation the description that uses 
atomic many-electron quantum numbers is most adequate; formally it is provided by the Green's function method for 
the Hubbard A-operators. 

In contrast with usual itinerant-electron ferromagnets, additional collective excitation branches (orbitons) occur. 
Due to the smallness of the current carrier concentrations, these modes possess low frequencies in comparison with 
typical electron energies such as the crystal-field splitting. Typically, the orbiton energies are of the same order 
of magnitude as the magnon energies. Also, mixed excitations with the simultaneous change of spin and orbital 
projections exist ("optical magnons"). All these excitations can be well defined in the whole Brillouin zone, the 
damping due to the interaction with current carriers being small enough. 

We have calculated one-particle Green's functions in the non-degenerate saturated ferromagnetic state. The ex- 
pressions obtained yield different criteria for spin and orbital instabilities. It turns out that the saturated spin 
ferromagnetism is more stable than the orbital one in the realistic case S > 1/2 (e.g., for magnetite and for colossal 
magnetoresistance manganites). This means that the half-metallic ferromagnetic phases both with saturated and 
non-saturated orbital moments can arise. A more detailed investigation with the use of real band and spectrum is 
required for concrete cases where orbital degrees of freedom yield important contributions to mag netic properties. In 
some situations (e.g., near the quantum phase transitions like those in virtual ferroelectrics 31]) the orbital modes 
can become soft, which should result in orbital instabilities at smaller conduction electron concentrations c. 

Due to the electron-orbiton interaction, additional contributions to the non-quasiparticle density of states can arise. 
Since the magnons are gapless, the electron-magnon contribution leads to a rather sharp one-sided increase of the 
density of states (with a crossover energy scale of the order of a typical magnon energy ZZJ) starting just from the 
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Fermi energy, Sg(E) oc (\Ep — E\/lJ) 3 ^ 2 [TH I23L |24| . In the case of half-metallic ferromagnets with unquenched 
orbital moment one can predict additional contributions to g(E) which correspond to orbital channels. Probably, the 
most simple way to probe the density of non-quasiparticle states is to measure the spin-polarized tunneling current 
[32t 133L l34j . Therefore it would be very interesting to investigate the density of states of "suspicious" materials by, 
e.g., spin-polarized scanning tunneling microscopy (STM) technique 35], the non-quasiparticle states being observable 
below the Fermi energy for the case of magnetite and above it for the case of colossal magnetoresistance manganites. 
It is worthwhile to mention that the non-quasiparticle states should exist also at the surface of the half-metallic 
ferromagnets 36] since STM is a surface-sensitive method. 
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(Support of Scientific Schools) , by the Russian Science Support Foundation and by the Netherlands Organization for 
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FIG. 1: Density of states in a half-metallic ferromagnet with double orbital degeneracy, I — —0.3 < for the semielliptic bare 
band with the width of W — 2, S = 1/2, crystal-field splitting A = 0.1. The energy is referred to the Fermi level. Solid and 
dashed lines correspond to crystal-field split subbands, the corresponding Fermi energies calculated from the band bottoms 
being 0.15 and 0.2, respectively. The spin-down subbands (lower half of the figure) nearly coincide with the bare bands shifted 
by IS. Non-quasiparticle states in the spin-up subbands (upper half of the figure) occur below the Fermi level. 
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FIG. 2: Density of states in a half-metallic ferromagnet with / = 0.3 > 0. The Fermi energies calculated from the orbital 
subband bottoms are 0.15 and 0.1. Other parameters are the same as in Fig.l The spin-up subbands (upper half of the figure) 
coincide with the bare bands shifted by —IS. The non-quasiparticle tails of the spin-down subbands (lower half of the figure) 
occur above the Fermi level. 
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FIG. 3: Density of states g,j m {E) in a half-metallic ferromagnet in the model l|13|l . The Fermi energy calculated from the 
bare band bottom is 0.2 (c = 0.034). Short-dashed line is the bare semielliptic density of states (a — +,m — +). The solid 
and dashed lines with peaks near the Fermi level show the contributions from orbital (a = +, m ^ +, A = 0.4) and spin 
(a = — ,ra = +) channels, respectively. The lower curves (without peaks) show the corresponding functions — (l/7r)ImR CTm (E). 
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FIG. 4: Density of states in a half-metallic ferromagnet in the model 11211 . The Fermi energy calculated from the bare band 
bottom is 0.1 (c = 0.019). Solid and short-dashed lines correspond to spin down electrons (the Green's function 1531 ) and spin 

up electrons (bare density of states). The dashed line shows the function — (l/7r)ImR (Tm (E). 



